Chiral effective field theory predictions for muon capture on deuteron and 3 He 
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The muon-capture reactions 2 H(p~, u^)nn and 3 He(/i _ , tv) 3 H are studied with nuclear potentials 
and charge-changing weak currents, derived in chiral effective field theory. The low-energy constants 
(LEC's) cd and ce, present in the three-nucleon potential and (cd) axial-vector current, are con- 
strained to reproduce the A — 3 binding energies and the triton Gamow- Teller matrix element. The 
vector weak current is related to the isovector component of the electromagnetic current via the 
conserved-vector-current constraint, and the two LEC's entering the contact terms in the latter are 
constrained to reproduce the A = 3 magnetic moments. The muon capture rates on deuteron and 
3 He are predicted to be 399 ± 3 sec -1 and 1494 ± 21 sec -1 , respectively. The spread accounts for 
the cutoff sensitivity as well as uncertainties in the LEC's and electroweak radiative corrections. By 
comparing the calculated and precisely measured rates on 3 He, a value for the induced pseudoscalar 
form factor is obtained in good agreement with the chiral perturbation theory prediction. 

PACS numbers: 23.40.-s,21.45.-v,27.10.+h 



When negative muons pass through matter, they can 
be captured into high-lying atomic orbitals. They then 
quickly cascade down into the IS orbit, where two com- 
peting processes occur: one is ordinary decay fj," — > 
e~ V e Vp, and the other is (weak) capture by the nucleus 
fi~ A(Z, N) Ujx A(Z - 1, N + 1). Apart from tiny cor- 
rections due to bound-state effects (chief among which is 
time-dilation) [H, the decay rate is essentially the same 
as for a free muon and, in light nuclei, is much larger than 
the rate for capture. The latter proceeds predominantly 
through the basic process p pT — )• n induced by the ex- 
change of a W + boson, and its rate, which would naively 
be expected to scale with the number Z of protons in the 
nucleus, is enhanced by an additional flux factor of Z 3 , 
originating from the square of the atomic wave function 
(w.f.) evaluated at the origin Q. Thus capture, with a 
rate proportional to Z 4 , dominates decay at large Z. 

Muon capture on hydrogen is, in principle, best suited 
to obtain information on the nucleon matrix element of 
the charge-changing quark current ^7^(1 — 75 )w, respon- 
sible for the process p pT — > nv^. Ignoring contributions 
from second-class currents Q for which there is presently 
no firm experimental evidence [H, it is parametrized in 
terms of four form factors (f.f.'s): two, F\{q 2 ) and i 7 ^? 2 ), 
from the polar-vector component of the weak current 
are related to the (isovector) electromagnetic (EM) form 
factors of the nucleon by the conserved-vector-current 
(CVC) constraint; two, the axial and induced pseu- 
doscalar f.f.'s G A (q 2 ) and Gps(q 2 ), go along with the 
axial- vector part of the weak current. The Fi(q 2 ) and 
F2(q 2 ) f.f.'s are well known over a wide range of momen- 
tum transfers q 2 from elastic electron scattering data on 
the nucleon @. Thef.f. G A (q 2 ) is also quite well known: 
its value at vanishing q 2 , g A = 1.2695(29), is from neu- 



tron /3-decay , while its g 2 -dependence is parametrized 



as G A (q 2 ) = g A /(l - q 2 /A 2 A ) 2 , with A A = 1 GeV from 
an analysis of pion electro-production data J7J and direct 
measurements of the reaction p — y n p + [8[ . 

Of the four f.f.'s, the induced pseudoscalar Gps(q 2 ) 
is the least known. The MuCap collaboration at PSI 
has recently reported a precise measurement of the rate 
for capture on hydrogen in the 15 singlet hyperfine 
state: r( 1 H)| singlct = 725.0 ± 13.7(stat) ± 10.7(syst) 
sec -1 Q. Based on this value, an indirect "experi- 
mental" determination of Gps at the momentum trans- 



fer q 2 = —0.88 m 2 relevant for p capture on hydrogen, 
Gpg P (<7q) = 7.3±1.2, has been obtained by using for the 
remaining f.f.'s the values discussed above and by eval- 
uating electroweak radiative corrections [l(J. The latter 
lead to a 2.8% increase in the rate on hydrogen, and are 
crucial for bringing Gpg P within less than la of the most 



recent theoretical prediction, G™(g 2 ) = 8.2 ± 0.2 [ill, 
obtained in chiral perturbation theory (%PT). For a re- 
cent and comprehensive review of theoretical and exper- 



13] 



imental efforts to determine Gps(q 2 ) see Refs. [1J 

In the present letter, we focus on the reactions 



2 R(p~ ,v^)nn and 3 He(//~, ^ M ) 3 H, hereafter referred to 
as p-2 and /i-3, respectively There are a couple of rea- 
sons for undertaking this study now: (i) the forthcoming 
measurement of the /x-2 rate r( 2 H) in the doublet hy- 
perfine state by the MuSun collaboration at PSI with a 
projected 1% precision [3, This and the already 
available, and remarkably precise, measurement of the 
fx-3 rate, T( 3 Hc) = 1496 ±4 sec" 1 0, will make it pos- 
sible to put tight constraints on Gp^ F (q 2 ) and to test the 
%PT prediction for this f.f. far more sharply than up to 
now. (ii) A number of low-energy weak processes of as- 
trophysical interest, such as the weak captures on proton 



2 



and 3 He, and neutrino reactions on light nuclei, are not 
accessible experimentally. In order to have some level 
of confidence in the reliability of their cross section esti- 
mates, it becomes crucial to study, within the same the- 
oretical framework, related electroweak reactions, whose 
rates are known experimentally, like muon captures flij ]. 
Theoretical work on the u-2 andu-3 reactions is quite 



extensive (see Refs. [12], |13|, LL8|)- So far, calcula 



tions have been performed within two different schemes: 
the "standard nuclear physics approach" (SNPA) and the 
approach known as "hybrid" chiral effective field theory 
(xEFT). In SNPA, Hamiltonians based on conventional 
two-nucleon (NN) and three-nucleon (NNN) potentials 
are used to calculate the nuclear w.f.'s, and the weak 
transition operator includes, beyond the one-body contri- 
bution (the impulse approximation — IA) associated with 



the basic process p [i 



meson-exchange currents 



as well as currents arising from the excitation of A-isobar 
degrees of freedom In the hybrid %EFT approach, 
the weak operators are derived in %EFT, but their ma- 
trix elements are evaluated between w.f.'s obtained from 
conventional potentials. Typically, the SNPA and hybrid 
%EFT predictions arc in good agreement with each other. 
For example, for the fi-2 rate, the SNPA calculation of 
Ref. [HI gives 391 sec -1 , to be compared with the hy- 
brid x EFT studies of Refs. [H and [3], which report 
386 sec -1 and 393 ± 1 sec" 1 , respectively. The differ- 
ences between Refs. [20[ and [HI are due to contributions 
of loop corrections and contact terms in the vector par t 
of the weak current, which were neglected in Ref. |20|. 
For the /x-3 rate, the SNPA calculation of Ref. [18j gives 
1486 sec" 1 , while the hybrid xEFT studies of Refs. [E| 
and [18( report, respectively, 1499±16sec _1 and 1484±4 



211 and hi 



sec . Here, the differences between Refs 
arise mostly from the inclusion in Ref. |21 of vacuum 
polarization effects on the muon bound state w.f. flp| — 
these would increase the SNPA and hybrid %EFT results 
of Ref. [l8| quoted above for the /x-3 rate to, respectively, 
1496 sec -1 and 1494 ± 4 sec" 1 . 

One of the objectives of the present work is to carry 
out a xEFT calculation of the fi-2 and /x-3 rates. Chi- 
ral EFT is a formulation of quantum chromodynamics 
(QCD) in terms of effective degrees of freedom suitable 
for low-energy nuclear physics: pions and nucleons. The 
symmetries of QCD, in particular its (spontaneously bro- 
ken) chiral symmetry, severely restrict the form of the in- 
teractions of nucleons and pions among themselves and 
with external electroweak fields, and make it possible 
to expand the Lagrangian describing these interactions 
in powers of Q/A x , where Q is pion momentum and 
A x ~ 700 MeV is the chiral-symmetry-breaking scale. 
As a consequence, classes of Lagrangians emerge, each 
characterized by a given power of Q / A x and each involv- 
ing a certain number of unknown coefficients, so called 
low-energy constants (LEC's). While these LEC's could 
in principle be determined by theory (for instance, in lat- 



tice QCD calculations), they arc in practice constrained 
by fits to experimental data. Some of them (for exam- 
ple, gA and the pion decay amplitude F v ) characterize 
the coupling (at lowest order) of pions to nucleons and, 
in particular, the strength of one- and two-pion-exchangc 
terms (denoted, respectively, OPE and TPE) in the NN 
potential (22l . [23j . i.e. its long-range components. Some of 
the other LEC's multiply NN (and multinucleon) contact 
interactions, and therefore encode short-range physics, 
which in a meson-exchange picture would, for example, 
be associated with vector-meson exchanges or excitation 
of baryon resonances, like the A isobar. 

The NN potential has been derived up to order 



(Q/A x ) 4 in the chiral expansion 22j, l23|. It consists 



of OPE and TPE with interaction vertices from lead- 
ing, next-to-leading, and next-to-next-to-leading 7rN chi- 
ral Lagrangians, and of contact terms. The LEC's have 
been constrained by accurate fits to the NN scattering 
database at e nerg ies below the pion production thresh- 
old (see Ref. [23[ for a review). The NNN potential, 
which first contributes at order (Q/A x ) 3 , includes S- and 
P-wave TPE — its P-wave piece is the familiar Fujita- 
Miyazawa NNN potential — a OPE plus NN contact term 
with LEC cd and a NNN contact terms with LEC ce- 

The vector and axial pieces of the weak current 
have been derived up to order Q/A x in, respectively, 



Refs. [2J, [25( and [26j. The one-body operators are the 



same as those obtained in the SNPA by retaining, in the 
expansion of the covariant single- nucleon four-current, 
corrections up to order (v/c) 2 relative to the leading- 
order term [191 ]. Two-body operators in the axial cur- 
rent (charge) first enter at order (Q/A x )° [(Q/A x ) -1 ], 
and are suppressed, in the power counting, by (Q/A x ) 3 
[Q/A x ] relative to the one-body term of order (Q/A x )~ 3 
[(Q/A x )~ 2 ]. In the axial current, these terms include a 
OPE contribution, involving the known LEC's C3 and C4 
(determined by fits to the NN data |23[), and one contact 
current, whose strength is parametrized by the LEC da 
(see below). In the axial charge, only OPE contributes, 
and the associated operator is proportional to gA/F 2 . 
One-loop corrections to the axial charge and current from 
TPE, which enter at Q/A x and are therefore strongly 
suppressed relative to the leading-order one-body terms, 
are ignored, since their contributions are expected to be 
tiny. 

The vector weak current is related (via the CVC con- 
straint) to the EM current, which includes, up to or- 
der Q/A x , OPE and TPE (i.e., one- loop corrections), 
as well as isoscalar and isovector contact terms, whose 
strengths are parametrized by the LEC's denoted, re- 
spectively, as g<is and g^v hi the following [24, 26 1. It has 
been shown [251 ] that such a current satisfies the continu- 
ity equation with the NN potential at order (Q/A x ) 2 . In 
this regard, we note that the construction of a conserved 
current with the (Q/A x ) 4 NN potential used here would 
require the inclusion of terms up to order (Q/A x ) 3 , i.e., 



3 



two-loop corrections. This is a daunting task, well be- 
yond the present state of the art. In a more speculative 
vein, it is also not obvious that such a theory could be 
made predictive, given the presumably large number of 
contact terms with unknown LEC's that it would entail. 

Finally, we notice that potentials and currents have 
power-law behavior for large momenta, and need to be 
regularized. This is accomplished in practice by introduc- 
ing a momentum-cutoff function. In the present work, 
the cutoff A is taken to be 500 MeV and 600 MeV. 

We now turn our attention to the determination of 
the LEC's dn, cd, ce, <?4s, and g^v- In the past, cd 
and ce were fixed by fitting the triton binding energy 
(BE) together with an additional strong-interaction ob- 
servable, such as the nd doublet scattering length 2 a n d 
or 4 He BE. However, this led to significant uncertain- 
ties, due to strong correlations between these observ- 
ables 27 1 . As the authors of Ref. [28| have observed, 
the LEC's dn and cd are related to each other via 
= T^fl CD + l M Jv(c 3 + 2c 4 ) + | (M N is the nu- 



cleon mass), and therefore one can fix cd (or dp) and 
ce by fitting the triton BE and half-life (specifically, 
the Gamow- Teller matrix element). Thus, we proceed 
as follows. The 3 H and 3 He ground state w.f.'s are cal- 
culated with the hyperspherical-harmonics method (see 
Ref. [H] for a review) using the chiral NN+NNN po- 



tentials of Rcfs. [22j, |23J, |30J for A = 500 and 600 MeV. 
The corresponding set of LEC's {cd,ce} is determined 
by fitting the A = 3 experimental BE's, BE( 3 H)=8.475 
MeV and BE( 3 He)=7.725 MeV, corrected for small con- 
tributions (+7 kcV in 3 H and -7 kcV in 3 He) due to 



the n-p mass difference [31[, since this effect is neglected 
in the present calculations. We then span the range 
Cd £ [— 3, 3], and, in correspondence to each cd in this 
range, determine ce so as to reproduce either BE( 3 H) or 
BE( 3 He). The resulting trajectories are shown in Fig. [1] 
and are nearly indistinguishable. Their average, shown 
by the red lines in Fig. [TJ leads to A = 3 BE's within 10 
keV of the experimental values above. Then, for each set 
of {cd, ce}, the triton and 3 He w.f.'s are calculated and, 
using the %EFT weak axial current discussed above, the 
Gamow- Teller matrix element of tritium /3-decay (GT TH ) 
is determined. The ratio GT TH /GT EXP is shown in 
Fig. [2j for both values of the cutoff A. We have used 

QrpEXP 



0.955 ± 0.004, as obtained in Ref. 
cept that we have conservatively doubled the error, rep- 
resented by the shadowed band in the figure. The range 
of c D values, for which GT TH = GT EXP within the ex- 
perimental error, arc [-0.20, -0.04] for A = 500 MeV, 
and [-0.32, -0.19] for A = 600 MeV. The corresponding 
ranges for c E are [-0.208,-0.184] and [-0.857,-0.833], 
respectively. We note that, for each pair of {cd,ce} in 
the selected range, the scattering length 2 a n d is calcu- 
lated to be 2 a nd = 0.666 ± 0.001 fm for A = 500 MeV 
and 2 a nd = 0.696 ± 0.001 fm for A = 600 MeV, which 
should be compared with 2 a nc i = 0.675 fm (29|, obtained 
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FIG. 1: (Color online) cd-ce trajectories fitted to reproduce 
the experimental 3 H and 3 He BE's. See text for explanation. 
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FIG. 2: The ratio GT TH /GT EXP as function of the LEC Co. 



with A = 500 MeV and {c D ,c E } = {1.0,-0.029}, as 
originally set in Ref. (131 . The most recent experimental 
determination gives 2 a n d = 0.645 ± 0.010 fm [321 ] . 

For the minimum and maximum values of {cd,ce} 
in the selected range, i.e., {—0.20,-0.208} and 
{-0.04, -0.184} for A = 500 MeV, and {-0.32, -0.857} 
and {-0.19, -0.833} for A = 600 MeV, we have deter- 
mined the isoscalar and isovector LEC's, g^s and gw, 
entering the NN contact terms of the EM current by re- 
producing the A = 3 magnetic moments. These LEC's 
are listed in Table HI 

TABLE I: The LEC's g4,s and g^v associated with the 
isoscalar and isovector NN contact terms in the EM current 
for A = 500 and 600 MeV. See text for explanation. 



{cd, ce} 



am 



8iv 



A=500 MeV {-0.20, -0.208} 0.207 ± 0.007 0.765 ± 0.004 
{-0.04,-0.184} 0.200 ±0.007 0.771 ± 0.004 

A=600 MeV {-0.32,-0.857} 0.146 ±0.008 0.585 ± 0.004 
{-0.19,-0.833} 0.145 ±0.008 0.590 ± 0.004 

Having fully constrained the NNN potential and weak 
current, we present in Table HJ the xEFT predictions for 
the fi~2 and fjr-S rates, T( 2 H) and T( 3 He). For T( 2 H), 
we also show the individual contributions of nn chan- 
nels with total angular momentum J < 2 ( 1 S'o, 3 Po, 3 -Pi, 
3 P2, 1 D2 and 3 i 7 2)- Higher partial waves are known to 
contribute less than 0.5 % to T( 2 H) jlj]. The one-body 
(IA) and (one+two)-body (FULL) results are listed sep- 



4 



arately. Note that the IA results depend on the cutoff 
A through the nuclear potentials. Theoretical errors in 
the FULL results arise from the fitting procedure, and are 
due primarily to the experimental error on 

GT EXP xhcy 

are not indicated when less than 0.1 sec -1 . Electroweak 
radiative corrections have been included as estimated in 
Ref. for hydrogen and 3 He — we have assumed that 
those for deuterium are the same as for hydrogen. The 
cutoff dependence of the predictions is weak, at less than 
1% level, thus suggesting that the mismatch between the 
chiral order of the potentials and that of the currents 
may be of little numerical import. If we also account 
for uncertainties in the electroweak radiative corrections 



of the order of 0.4% [l0|, we can conservatively quote 
r( 2 H) = 399 ± 3 sec" 1 and r( 3 He) = 1494 ± 21 sec" 1 . 
These predictions are in good agreement with available 
experimental data (although those on r( 2 H) [h| have 
large errors), as well as with results of recent theoreti- 
cal studies pil [2(| [2l|. Finally, a comparison between 
the calculated and measured n~3 rates makes it possi- 
ble to put a constraint on the induced pseudoscalar f.f. 
Gps(q 2 ) at g 2 = —0.954m 2 relevant for the /i-3 reac- 
tion. By varying Gps{q$) so as to match the theoretical 
upper (lower) value with the experimental lower (upper) 
value for the rate, we obtain Gps(g 2 ) = 8.2±0.7, in good 
agreement with the xPT prediction of 7.99 ± 0.20 [111 ]. 



TABLE II: Total rates for muon capture on deuteron F( 2 H) and 3 He F( 3 He), in sec -1 , corresponding to A = 500 and 600 MeV. 
The one-body (IA) and (one+two)-body (FULL) contributions are listed, along with the individual partial-wave contributions 
to r( 2 H). Theoretical uncertainties in the FULL results, not reported when below 0.1 sec -1 , are due to the fitting procedure. 





'So 6 Po ''Pi A P2 'D 2 *P 2 T( 2 H) 


r( 3 He) 


IA(A = 500 MeV) 
IA(A = 600 MeV) 


238.8 21.1 44.0 72.4 4.5 0.9 381.7 
238.7 20.9 43.8 72.0 4.5 0.9 380.8 


1362 
1360 


FULL (A = 500 MeV) 
FULL (A = 600 MeV) 


254.4T0.9 20.5 46.8 72.1 4.5 0.9 399.2T0.9 
255.2T1.0 20.3 46.6 71.6 4.5 0.9 399.1T1.0 


1488±9 
1499±9 
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